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Abstract 

Using constituent quark model constraints we calculate the gluon and sea-quark 
content of pions solely in terms of their valence density (fixed by 7riV Drell-Yan data) 
and the known sea and gluon distributions of the nucleon, using the most recent 
updated valence-like input parton densities of the nucleon. The resulting small-x 
dynamical QCD predictions for g n (x,Q 2 ) and q n (x,Q 2 ) are unique and parameter 
free. Simple analytic parametrizations of the resulting parton distributions of the 
pion are presented in LO and NLO. These results and parametrizations will be 
important, among other things, for updated formulations of the parton distributions 
of real and virtual photons. 



The parton content of the pion is poorly known at present. The main experimental 
source about these distributions is mainly due to data of Drell-Yan dilepton production 
in 7r~-tungsten reactions ^ |3[], which determine the shape of the pionic valence density 
v 7T (x, Q 2 ) rather well, and due to measurements of direct photon production in n ± p 
7X [I], |J which constrain the pionic gluon distribution g' K {x 1 Q 2 ) only in the large-a; 
region ||. In general, however, present data are not sufficient for fixing g 77 uniquely, in 
particular the pionic sea density q n (x, Q 2 ) remains entirely unconstrained experimentally. 
Therefore we have previously || utilized a constituent quark model to relate q n and 
g n to the much better known radiatively generated parton distributions f p (x, Q 2 ) of the 
proton ||. These relations arise as follows: describing the constituent quark structure 
of the proton p = UUD and the pion, say n + = UD, by the scale (Q 2 ) independent 
distributions U p,7T+ (x), D p (x) and Z) 7r+ (x), and their universal (i.e. hadron independent) 
partonic content by v c (x, Q 2 ), g c (x,Q 2 ) and q c (x,Q 2 ), the usual parton content of the 
proton and the pion is then given by 



f p (x,Q 2 ) = f 1 ^y[U p (y) + D p (y)]f c (-,Q 2 ) (1) 

J' y \y J 

F(x,Q 2 ) = l ld ^[u- + {y) + D^{y)]f c ( X -,Q 2 ) (2) 
Jx y 1 J \y J 

where f = v, q, g with v p = u p + d p , q p = (u p + d p )/2 J v« = < + + df, q n = 
( , u 7r+ + d 7T+ )/2 and w 7r+ = d w+ due to ignoring minor SU(2)fl avor breaking effects in the 
pion 'sea' distributions. Assuming these relations to apply at the low resolution scale 
Q 2 = (a*lo = 0.23 GeV 2 , 

A'nlo = 0.34 GeV 2 ) of || where the strange quark content 
was considered to be negligible, 

s p {x, fi 2 ) = s p {x, ft 2 ) = s ff (x, fi 2 ) = s*(x, fi 2 ) = 0, (3) 

one obtains from (1) and (2) the constituent quark independent relations 

v^jn,/! 2 ) = q w (n,fi 2 ) = g w (n,fi 2 ) 
v p (n,[i 2 ) q p (n,fi 2 ) g p (n,fi 2 ) 

where for convenience we have taken the Mellin n-moments of eqs. (1) and (2), i.e. 
f(n,Q 2 ) = Slx n ~ x f{x,Q 2 )dx. Thus, as soon as v n (x,fi 2 ) is reasonably well determined 



from experiment, our basic relations (4) uniquely fix the gluon and sea densities of the 
pion in terms of the rather well known parton distributions of the proton: 

9{n,H) = — p n,/i , q ( n , n ) = — -q p {n,(i). 5 



Furthermore, the sum rules 



l 

v w {x,Q 2 )dx = 2 (6) 

o 

xv^[x,Q 2 )dx = [ xv p (x,Q 2 )dx (7) 
10 Jo 

impose strong constraints on f ?r (x, fi 2 ) which are very useful for its almost unambiguous 
determination from the ttN Drell-Yan data. Notice that eq. (7), together with (4), implies 
the energy-momentum sum rule for f 7 * to be manifestly satisfied. In addition, eq. (7) 
implies that the valence quarks in the proton and the pion carry similar total fractional 
momentum as suggested by independent analyses within the framework of the radiative 
parton model pi RJ. 

The relations in eq. (5) imply that any updating of f p (x,fi 2 ) yields a corresponding 
updating of / 7r (x,/i 2 ). Recently an updating of / p (x,/x 2 ) within the framework of the 
radiative (dynamical) parton model was undertaken utilizing additional improved data 
on F p (x, Q 2 ) from HERA [|T0l , [TT|1 and a somewhat increased a s (M|) = 0.114 resulting in 



a slight increase in /i 2 (yU 2 Q = 0.26 GeV 2 , /i^Lo = 0-40 GeV 2 ). An improved treatment 
of the running a s (Q 2 ) at low Q 2 was furthermore implemented by solving in NLO(MS) 



numerically |9|] rather than using the approximate NLO solution 

a s (Q 2 ) 1 /3i lnln{Q 2 /A 



(9) 



4vr p ln(Q 2 /A 2 ) /5 3 In 2 (Q 2 /A 2 ) 

as done in || || H, which is sufficiently accurate only for Q 2 £ m 2 ~ 2 GeV 2 0. The 
LO and NLO evolutions of f n (n, Q 2 ) to Q 2 > /i 2 are performed in Mellin n-moment 



space, followed by a straightforward numerical Mellin-inversion [|12| to Bjorken-x space 



It should be noted that the evolutions are always performed in the fixed (light) / = 3 
flavor factorization scheme |L3|, |6|, || [J, i.e. we refrain from generating radiatively massless 
'heavy' quark densities h n (x,Q 2 ) where h = c, b, etc., in contrast to ||. Hence heavy 
quark contributions have to be calculated in fixed-order perturbation theory via, e.g., 
g-KgP fofo^ u^vP — > hh, etc. (Nevertheless, rough estimates of 'heavy' quark effects, 
valid to within a factor of 2, say, can be easier obtained with the help of the massless 
densities c^^x^Q 2 ) and ^(x, Q 2 ) given in ||.) 

Using all these modified ingredients together with the new updated |J f p (x,fi 2 ) in 
our basic predictions in eq. (5), the present reanalysis of the available Drell-Yan data 0, 
closely following the procedure described in 0, yields 

Ulo(z^lo) = l-129x- a496 (l-x) a349 (l + 0.153 v / ^) (10) 

^NLO(^NLO) = 1.391X- 0447 (1 -X) - 426 (11) 

where /Ilo = 0-26 GeV 2 and /i^Lo = 0-40 GeV 2 . These updated input valence densities 
correspond to total momentum fractions 

/ xu lo( x ; ^lo)^ 3 - = 0.563 (12) 
J o 

/ ^NLoO^NLC))^ = °- 559 ( 13 ) 

as dictated by the valence densities of the proton via eq. (7). Our new updated input 
distributions in eqs. (10), (11) and (5) are rather different than the original GHV^ input 
0] in fig. 1 which is mainly due to the vanishing sea input of GRV^ in contrast to the 
present one in eq. (5). On the other hand, our updated input in fig. 1 is, as expected, 
rather similar to the one of . In both cases, however, the valence and gluon distributions 
become practically indistinguishable from our present updated ones at scales relevant for 
present Drell-Yan dimuon and direct-7 production data, Q 2 = M 2 +fi - ~ 20 GeV 2 , as 
illustrated in fig. 2. Therefore our present updated pionic distributions give an equally 
good description of all available ttN Drell-Yan data as the ones shown in ||. 

For completeness let us mention that our basic predictions (5) for the valence-like 
gluon and sea densities at Q 2 = fi 2 , as shown in fig. 1, can be simply parametrized in 
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Bjorken-x space : in LO at Q 2 = fi1 Q = 0.26 GeV 2 

xg*{x,nlo) = 7.326 x L433 (l - 1.919 y/x + 1.524 x)(l-x) L326 
xq n (x,nl ) = 0.522 x ai60 (l - 3.243 ^x + 5.206 x)(l - x) 5 ' 20 , (14) 

whereas in NLO at Q 2 = /!n LO = 0.40 GeV 2 we get 

xg' K {x 1 ^ Q ) = 5.90 x L270 (l- 2.074 + 1.824x)(l-x) L290 
xq^fafj&uo) = 0.417 x a207 (l - 2.466 jx + 3.855 x) (1 - x) 4 ' 454 . (15) 



Finally, fig. 3 shows our resulting predictions for xg n (x,Q 2 ) and xq n (x,Q 2 ) as com- 
pared to the former GRV,,- results [Q. The GRV„- results for xq 7 " are significantly steeper 
and softer for x £ 0.01 due to the vanishing SU(3)fl avor symmetric (light) sea input 
xg"(x, l u 2 ) = 0, in contrast to our present approach || based on a more realistic fi- 
nite light sea input in eq. (5). The valence-like gluon and sea inputs at Q 2 = fi 2 , which 
become (vanishingly) small at x < 10~ 2 , are also shown in fig. 3. This illustrates again 
the purely dynamical origin of the small-x structure of gluon and sea quark densities at 
Q 2 > fi 2 . Our predictions for s n = s n , as evolved from the vanishing input in eq. (3), are 
not shown in the figure since they practically coincide with q n (x, Q 2 ) of GRV,,- shown in 
fig. 3 which also results from a vanishing input ||. Simple analytic parametrizations of 
our LO and NLO predictions for f n (x, Q 2 ) are given in the Appendix. 

To conclude let us recall that an improvement of f^iXjQ 2 ) is particularly important 
in view of its central role in the construction of the photon structure function and the 



photonic parton distributions |L4|, [15], [16], [L7|, [Tj|. Furthermore, recent (large rapidity gap) 



measurements of leading proton and neutron production in deep inelastic scattering at 
HERA [IS] allow, under certain (diffractive) model assumptions, to constrain and test the 



pion structure functions for the first time at far smaller vales of x (down to about 10 3 ) 
than those attained from fixed target ttN experiments. 
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Appendix 

A. Parametrization of LO parton distributions 
Defining M 



In 



In [Q7(0.204 GeV) 2 ; 
In [/i 2 o /(0.204 GeV) 2 



(A.l) 



to be evaluated for /j^ = 0.26 GeV 2 , all our resulting pionic parton distributions can 
be expressed by the following simple parametrizations, valid for 0.5 ^ Q 2 ;$ 10 5 GeV 2 
(i.e. 0.31 < s < 2.2) and 10~ 5 5, x < 1. For the valence distribution we take 



xv 



w {x, Q 2 ) = Nx a {l + + Bx)(l - x) D 



(A.2) 



with 



N 
a 
A 
B 
D 



1.212 + 0.498 s + 0.009 s 2 
0.517 -0.020 s 
-0.037 - 0.578 s 
0.241 + 0.251s 
0.383 + 0.624 s. 



(A.3) 



The gluon and light sea-quark distributions are parametrized as 



x a (A + B^ + Cx) (in ~) + s a exp (-E + ^E'sHn^ 



1 — x 



(A.4) 
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For w = g 



a = 0.504, $ = 0.226 



a 



= 2.251- 1.339 Vs, 6=0, 

A = 2.668- 1.265 s + 0.156 s 2 , B = -1.839 + 0.386 s, (A.5) 

C = -1.014 + 0.920 s -0.101s 2 , D = -0.077 + 1.466 s, 

E = 1.245 + 1.833 s, E' = 0.510 + 3.844 s , 



and for the light sea w = q 



a 


= 1.147, 




= 1.241, 


a 


= 0.309- 0.134 v/i, 


b 


= 0.893 - 0.264^ 


A 


= 0.219 - 0.054 s, 


B 


= -0.593 + 0.240 s. 


C 


= 1.100 - 0.452 s, 


D 


= 3.526 + 0.491s, 


E 


= 4.521 + 1.583 s, 


E' 


= 3.102. 



The strange sea distribution s 71 " = s 77 is parametrized as 



— 7T / 



a = 0.823, (3 = 0.650, 

a = 1.036 -0.709 s, A = -1.245 + 0.713 s, 

B = 5.580 - 1.281s, D = 2.746 - 0.191 s, 

E = 5.101 + 1.294 s, E' = 4.854 - 0.437 s . 



(A.6) 



xs " [x, Q 2 ) = jj~Ty (l + + Bx)(l- x) D exp ^—E + y E'sHn^ J (A.7) 

with 



(A.8) 



B. Parametrization of NLO(MS) parton distributions 
Defining 

Zn[Q 2 /(0.299GeV) 2 ] 
S " m in [^ LO /(0.299 GeV) 2 ] ^ 

to be evaluated for ^ LO = 0.40 GeV 2 , our NLO predictions can be parametrized as the 
LO ones and are similarly valid for 0.5 <, Q 2 < 10 5 GeV 2 (i.e. 0.14 < s & 2.38) and 
10~ 5 £ x < 1. The valence distribution is given by (A. 2) with 

N = 1.500 + 0.525 s -0.050 s 2 
a = 0.560 - 0.034 s 
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A = -0.357 - 0.458 s 
B = 0.427 + 0.220 s 

D = 0.475 + 0.550 s. (A.10) 

The gluon and light sea distributions are parametrized as in (A. 4) where for w = g 

a = 0.793, = 1.722, 

a = 1.418-0.215^5, 6=0, 

A = 5.392 + 0.553 s - 0.385 s 2 , B = -11.928 + 1.844 s, (A.ll) 

C = 11.548 - 4.316 s + 0.382 s 2 , D = 1.347 + 1.135 s, 

E = 0.104+ 1.980 s, E' = 2.375 - 0.188 s . 

and for the light sea w = q 

a = 1.118, = 0.457, 

a = 0.111 - 0.326 Vs, b = -0.978 - 0.488^, 

A = 1.035 - 0.295 s, B = -3.008 + 1.165 s, (A.12) 

C = 4.111 - 1.575 s, D = 6.192 + 0.705 s, 

E = 5.035 + 0.997 s, E' = 1.486 + 1.288 s . 

The strange sea distribution is parametrized as in (A. 7) with 

a = 0.908, p = 0.812, 

a = -0.567 - 0.466 s, A = -2.348 + 1.433 s, 

B = 4.403, D = 2.061, 

E = 3.796 + 1.618 s, E' = 0.309 + 0.355 s . 



(A.13) 



Let us recall that in the light quark sector u*J = = = , u n+ = d n+ = 
u *- = d*~ and p° = (/ 7r+ + p~)/2. 
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Figure Captions 



Fig. 1 The valence and valence-like input distributions xf 7T (x, Q 2 = /j, 2 ) with / = v, q, g as 
compared to those of GRVV ||. Notice that GRV^ employs a vanishing SU(3)fl avor 
symmetric g 71 " input at yU 2 Q = 0.25 GeV 2 and /Un LO = 0-3 GeV 2 p. Our present 
SU(3)fl avor broken sea densities refer to a vanishing s n input in (3), as for GRV,,- ||. 

Fig. 2 Comparison of our NLO valence distribution at Q 2 = 20 GeV 2 with the one of 
GRV^ and GRS ||. This density plays the dominant role for describing presently 
available nN Drell-Yan dimuon production data. For illustration, the gluon and 
sea densities are shown as well. The SU(3)fl avor symmetric GRV,,- sea q n = s n is 
not shown, since it is similar to s n of our present analysis and of GRS which are all 
generated from a vanishing input at Q 2 = /i 2 , cf. eq. (3). 

Fig. 3 The small-x predictions of our radiatively generated pionic gluon and sea-quark 
distributions in LO and NLO at various fixed values of Q 2 as compared to those of 
GRV,r 0. The valence-like inputs, according to eq. (5) as presented in fig. 1, are 
shown for illustration by the lowest curves referring to fi 2 . The predictions for the 
strange sea density s 71 = s 71 are similar to the GRV,,- results for q n . The results are 
multiplied by the numbers indicated in brackets. 
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Fig. 1 




Fig. 2 



